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We present a systematic study on the excited electron-bubble states in superfluid 4 He using a 
time-dependent density functional approach. For the evolution of the IP bubble state, two different 
functionals accompanied with two different time-development schemes are used, namely an accurate 
finite-range functional for helium with an adiabatic approximation for electron versus an efficient 
zero-range functional for helium with a real-time evolution for electron. We make a detailed com- 
parison between the quantitative results obtained from the two methods, which allows us to employ 
with confidence the optimal method for suitable problems. Based on this knowledge, we use the 
finite-range functional to calculate the time-resolved absorption spectrum of the IP bubble, which 
in principle can be experimentally determined, and we use the zero-range functional to real-time 
evolve the 2P bubble for several hundreds of picoseconds, which is theoretically interesting due to 
the break down of adiabaticity for this state. Our results discard the physical realization of relaxed, 
metastable 2P electron-bubbles. 
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I. INTRODUCTION 

Electron bubble (e-bubble) in liquid-helium has been 
an attractive topic for numerous experimental and theo- 
retical studies in the past, and it has drawn again some 
interest in recent years GHH Density functional (DF) the- 
ory has proved to be a powerful tool in dealing with many 
interesting physical situations involving electron bubbles. 
When it is applied to optically excited e-bubbles, not only 
can it achieve quantitativ e agreem ent with experiments 
on the absorption spectra, ^ * 12 l 13 l it can also nicely dis- 
play the dynamical evolutions on the picosecond time 
scale, such as how the bubbles change shapes, release 
energy, or even break into smaller bubbles!- 14 * 15 * These 
latest works likely require using different density func- 
tionals (finite-range or zero-range for liquid-helium) in 
different time-development schemes (adiabatic or real- 
time evolutions for electron). Regardless of the technical 
details, the time-dependent density functional approach 
is no doubt the only workable approach at present for 
studying the evolution of the excited e-bubble states in 
liquid-helium. The quantitative results drawn from these 
simulations can be useful to interpret the experimental 
results and predict new ones. 

Upon dipole excitation from the IS ground state to 
the IP or 2P excited state, an e-bubble evolves by re- 
laxing its shape around the excited electron probability 
density. This relaxation eventually drives the e-bubble 
back to the spherical IS ground state. It has been quan- 
titatively showrP^ that depending on the pressure (P) 
applied to the liquid, this may happen in two different 
ways. At a pressure below about 1 bar, the e-bubble 
undergoes damped oscillations for a period of time long 
enough to allow the electron to radiatively decay to the 
deformed IS state, which then evolves radiationlessly to 



the spherical IS state. In contrast, at 1 bar or above, the 
excited e-bubble evolves towards a configuration made 
of two baby bubbles, so that the probability of finding 
the electron evenly distributes between them. This two- 
bubble configuration is unstable against asymmetric per- 
turbations, and one expects this instability to cause the 
electron to loc alize in one of the baby bubbles while the 
other collapses P 2 1 15 ' 16 1 

When doing this calculation, the authors of Ref. [15] 
naturally chose the so-called Orsay-Trento DFP^I This 
functional is finite-range and incorporates a term that 
mimicks back-flow effects in order to accurately repro- 
duce the dispersion relation of the elementary excitations 
of superfluid 4 He. This is instrumental to properly de- 
scribe the energy transfer from the bubble to the liquid, 
that proceeds by causing all sorts of possible excitations 
in the superfluid. The token one has to pay for its use is 
the very high computational cost. Due to the large dif- 
ference between the intrinsic time scales of electron and 
helium, this functional is not well adapted for fully real- 
time, three-dimensional evolutions. Considering this lim- 
itation, the adiabatic approximation was used to update 
the electron wavefunction at every instantaneous helium 
configuration.^ One of the main concerns is then to es- 
tablish how long the adiabatic approximation is valid for. 
A careful analysis led to the conclusion that it holds for 
the IP e-bubble for at least several tens of picoseconds 
a period of time large enough to guarantee the reliabil- 
ity of the results obtained for this state. Contrarily, the 
adiabatic approximation breaks down very quickly for 
the 2P e-bubble, implying that the existence of relaxed 
quasi-equilibrium 2P bubbles is questionable. 

The authors of Ref. [14] followed a different path. In- 
stead of using a finite-range DF, they employed a much 
simpler zero-range one. This allowed them to carry out 
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fully real-time calculations for hundreds of picoseconds 
without imposing any adiabatic assumption, at the price 
of an inaccurate description of the elementary excitations 
of the liquid. Although this may not qualitatively affect 
the physical results, it is unclear how quantitatively reli- 
able are the results so obtained. 

In this work, we perform a more systematic study on 
the relaxation of the excited e-bubbles. We first carry 
out a detailed comparison between the finite-range and 
zero-range functionals applied to the IP bubble problem, 
through which we gain some insight about their strong- 
points and shortcomings. We then use the finite-range 
functional to calculate the time-resolved absorption spec- 
trum of the IP bubble, which can in principle be mea- 
sured in the experiments. Next, we use the zero-range 
functional to real-time evolve the 2P bubble for several 
hundreds of picoseconds, which is of theoretical interest 
as one can clearly trace how the adiabatic approximation 
breaks down for this state. 

Our paper is organized as follows. In Sec. II, we intro- 
duce the theoretical framework and numerical schemes. 
In Sec. Ill, we discuss the results so obtained. Finally, a 
summary is presented in Sec. IV. Several movies showing 
the dynamical evolution of electron bubbles can be found 
in the supplementary material.^ 



II. THEORETICAL FRAMEWORK 

Within the DF approach, the energy of an electron- 
helium system at zero temperature can be written as a 
functional of the single-electron wavefunction <I> and the 
macroscopic helium wavefunction ^ 
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Specifically, ^ = ^pexp[zS] gives the helium particle 
density p and the superfTuid velocity v = WS/mH e . 
V^_He[p] i s the electron- helium interaction potential,^ 
and f He-He [p] is the helium- helium potential energy 
density. For the sake of comparison, we choose for 
^He-He[p] either the finite-range Orsay-Trento (OT) den- 
sity functionalp'or the zero-range Stringari-Treiner (ST) 
density functional. 20 It is known that the former one pro- 
vides a very accurate description of superfluid 4 He, par- 
ticularly of the dispersion relation that covers all the ro- 
ton excitations up to the wave number q = 2.3 A -1 . 
In contrast, the latter one only reproduces the long- 
wavelength phonon excitations and is not so accurate, 
but it has the advantage of high computational efficiency 
in dynamic evolutions. The dispersion relations obtained 
from both functionals using the method of Ref. [21] are 
plotted in Fig. [I] together with the experimental results.^ 
Functional variations of the associated grand poten- 
tial with respect to \I/ and $ yield the following Euler- 



Lagrange equation for the helium and Schrodinger equa- 
tion for the electron: 
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where \i is the helium chemical potential and e is the 
electron eigenenergy. Throughout this paper, we treat 
pressure as a given external condition. The associated 
chemical potential and particle density in bulk liquid are 
obtained from the equation of state derived from the DF 
being used. 

The above equations are solved numerically with 13- 
point finite-difference formulae. In the finite-range DF 
case, we work in three-dimensional Cartesian coordinates 
that allow for an extensive use of fast Fourier transfor- 
mation techniques 23 as explained in Ref. El Whenever 
necessary, we implement a Gram- Schmidt scheme to de- 
termine from Eq. Q the electron spectrum in the helium 
cavity. In the zero-range DF case, we work in cylindrical 
coordinates assuming azimuthal symmetry around the z- 
axis (r = 0) and specular symmetry about the z = 
plane. Hence, we only need to solve the equations in the 
r > and z > "quadrant", which greatly speeds up 
the calculation and so we can attain longer-time physics. 
There is no difficulty to relax these symmetry restric- 
tions except slowing down the calculation. However, the 
results are not expected to be very different since there is 
no detectable symmetry-breaking instability in our prob- 
lems. For both functionals we use a fairly large spatial 
step of about 1 A, without any apparent loss of numerical 
accuracy.^ 

The time evolution starts from an excited e-bubble 
state, which means that the electron has been suddenly 
brought from the IS onto the IP or 2P state of the orig- 
inal spherical e-bubble. From this initial configuration, 
the superfluid helium then evolves according to 
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where the detailed form of the effective potential 
^He-He[p, v ] can be found e.g. in Refs. [24landl25l 

In the adiabatic approximation scheme, we do not 
evolve the electron in real-time but keep tracing the in- 
stantaneous eigenstates satisfying Eq. (J3|. In most cases, 
e-bubbles around excited electron states evolve towards 
configurations that are not spherically but axially sym- 
metric. The originally degenerate angular momentum 
electron eigenstates in the spherical bubble now split ac- 
cording to the projection (m-values) on the symmetry 
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z-axis, among which only the ±ra states are still de- 
generate with each other. It is thus convenient to use 
the notation for the orbital angular momentum of sin- 
gle particle states in linear molecules, i.e., cr, 7r, 5, 0, . . . 
for \m\ = 0,1,2,3,... In addition to the axial symme- 
try, a nP bubble also keeps the original specular symme- 
try in the course of its evolution. Hence, one can con- 
struct the electron eigenbasis in such a way that the elec- 
tron wavefunctions satisfy <I>(r, z, 0) = ±$(r, —z, 0). The 
correspondence between the lower lying spherically and 
axially symmetric electron states is displayed in Fig. [2] 
along with a representation of their probability densities. 
The superscript + (— ) denotes specularly symmetric (an- 
tisymmetric) states. The adiabatic evolution is obtained 
by keeping the IP (2P) electron in the instantaneous 1<j~ 
(3<j~) eigenstate. 

In the real-time dynamics scheme, the electron evolves 
according to 
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We employ a fourth-order Runge-Kutta method to ob- 
tain the first time steps, and Hamming's method^ for 
subsequent steps. A time step of 10 -2 ps is chosen for 
the adiabatic evolution, and of 10 -6 ps for the dynamical 
evolution. These very different values reflect the mass ra- 
tio m e /mH e ~ 10 -4 . This is this reason that makes a dy- 
namical evolution unaffordable when f He-He [p] is finite- 
range, since updating the mean field is computationally 
very costly. 

During the bubble evolution, sound waves released 
from its surface eventually reach the cell boundary. If 
no action is taken, they will bounce back spoiling the 
calculation. A way to handle this problem is to include 
some source of damping into Eq. Q governing the fluid 
evolution, see e.g. Refs. 1 27 28, We have opted by mak- 
ing the replacement i — > i + A(r) in Eq. Q. This 
corresponds to a rotation of time axis in the complex 
plane by introducing a damping field A(r), which takes 
the forrrP^ 
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We keep the dimensionless parameter Aq — 1.6, and set 
a = 5 A, so = 60 A in the finite-range calculation, and 
a = 8 A, sq = 90 A in the zero-range calculation. The 
evolution is damping-free [A(r) < 1] in a sphere of ra- 
dius s < so — 2a, which is 50 A in the finite-range case 
and 70 A in the zero-range case. From Figs. [5][6] and the 
supplementary material 18 one can see that this region is 
large enough for the IP e-bubble to expand within an 
undampening environment. For the 2P e-bubble evolu- 
tion, we use a (r, z) calculation box of 150 x 150 A 2 and 
so = 120 A, leaving ~ 100 A of undampening space for 
the bubble to expand. 

The above prescription works extremely well, as it effi- 
ciently dampens the excitations of the macroscopic wave- 



function at the cell boundaries, and does not need a large 
buffer region to absorb the waves — actually we use the 
same box where the starting static calculations have been 
carried out It allows us to extend the adiabatic calcu- 
lations of Ref. [15] from tens to hundreds of picoseconds. 



III. RESULTS AND DISCUSSION 
A. IP e-bubble dynamics 

1. Adiabatic versus real-time dynamical evolution 

To some extent, an e-bubble in liquid-helium is nearly 
a textbook example of an electron confined in a spherical- 
square-well potential. Its static properties are fairly in- 
sensitive to the complexities of the chosen functional 
provided t he bu lk and surface properties can be well 
reproduced! 12 * 13 * In particular, a zero-range DF descrip- 
tion of the IS- IP absorption energies of the e-bubble as a 
function of P does not differ much from that obtained by 
a finite-range DF description! 3 -^ This means that in some 
situations one may simply use a zero-range functional, 
which has an advantage of high computational speed. 

Clearly, the dynamics of an e-bubble is much more 
involved than its statics. In our problems, the - 
nonspherical- squeezing and stretching of the bubble may 
cause its waist to shrink to a point when electron tunnel- 
ing plays a role, and may also dissipate a large amount of 
energy by exciting elementary modes in the surrounding 
liquid. So, even if the static properties of the e-bubble 
are equally well described by both functionals, it is not 
obvious whether they yield a similar dynamical evolution 
for an excited e-bubble. This is the first issue we want to 
address. 

We use two different schemes to compute the relaxation 
of a IP e-bubble at P = 0, 0.5, 1 and 5 bars, and compare 
the results so obtained. One such scheme is the finite- 
range OT density functional description for helium with 
the adiabatic evolution for electron. The other scheme 
is the zero-range ST density functional description for 
helium with real-time evolution for the electron. 

As can be seen in Figs. [3][6j the evolution starts with 
the bubble stretching along the z-axis and shrinking on 
its waist. After this stage, the bubble may continue os- 
cillating and releasing energy into the liquid, eventually 
reaching a relaxed, metastable IP state, or may split into 
two baby bubbles due to the liquid filling-in around the 
bubble waist. The density waves radiated to the liquid 
during this evolution take away a considerable part of 
the energy injected into the system during the absorp- 
tionprocess, i.e., 105 meV at P = and 148 meV at 5 
barP 

In Figs. [3][6]we compare the bubble evolution obtained 
within the two frameworks for different pressures. At a 
first glance, both dynamics are nearly equivalent dur- 
ing the first 50 ps, starting to differ from this time on 
although they are still qualitatively similar. A more de- 
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tailed analysis, focused on three key elements of the den- 
sity profiles, indicates the following: 

a. The shape of the bubble surface, defined as the 
locus where the liquid density equals half the saturation 
density value p , e.g., 0.0218 A" 3 at P = bar. This 
shape determines the most crucial properties of an e- 
bubble. From this shape, we know whether the bubble is 
simply connected or has split. Up to t < 50 ps, the shape 
of the bubble is nearly identical in both descriptions. At 
later times, the bubble shape changes at a slower pace in 
the ST than in the OT description. 

Figure [7] illustrates the time evolution of the bubble 
surface. In particular, the top panel shows the evolu- 
tion of the point on the bubble surface at r = with 
z > 0. This represents half the longitudinal extent of the 
e-bubble. One can see that this length oscillates in the ST 
calculation with a lower frequency than in the OT one. 
If the bubble symmetrically splits into two baby bubbles, 
there are two such points, as can be seen in the bottom 
panel for P = 1 bar. We have found that, besides the 
moment at which the distance between the baby bubbles 
increases steadily — about 175 ps for the OT functional 
and 200 ps for ST functional — there is a time interval 
between ~ 60 and ~ 90 ps for the OT functional, and ~ 
60 and ~ 110 for the ST functional, where the IP bubble 
at P = 1 bar has split but the emerging baby bubbles 
are "locked" by the shared electron that exerts some at- 
tractive force on them, forcing them back to a simply 
connected configuration. Eventually, the baby bubbles 
are unlocked and the distance between them grows. 

b. The surface thickness of the bubble, defined as the 
width of the region satisfying O.lpo < p(r) < 0.9po- The 
thickness of the bubble surface has been found to be 
nearly independent of the local surface curvature at any- 
time during the evolution, see also Ref. [T5l It i s abo ut 
1 A lar ger in the ST than in the OT description) 11 as 
can be seen in Figs. [3][6] (the blurrier the bubble-helium 
interface, the larger the surface thickness). The zero tem- 
perature OT result, a bout 6 A, is in agreement with the 
experimental findings! 31 * 32 * 

c. The density oscillations traveling through the liq- 
uid. This is the point at which the differences between 
the two functionals become more apparent. The density 
waves produced by the ST functional have much larger 
wavelengths because this approach cannot sustain short 
wavelength inhomogeneities due to the huge energy cost 
fromthe |Vp(r)| 2 surface energy term. The OT func- 
tional has not such a term and is free from this drawback. 
Roughly speaking, the short wavelength waves arising in 
the OT approach are smeared out in a sort of big tsunami 
in the ST case, see for instance the panels corresponding 
to t = 5 and 10 ps in Fig. [6] It is worth emphasiz- 
ing that the wave interference pattern found in the OT 
description is not an artifact produced by waves bounc- 
ing back from the box boundaries, as those are already 
washed out by the damping term. It arises from the in- 
terference of waves produced at different points on the 
e-bubble surface. 



To quantitatively study the nature of the waves emit- 
ted during the bubble evolution, we perform a Fourier 
analysis of the density profile along the symmetry axis, 
restricting it to the region 30 A < z < 70 A, away from 
the bubble location to avoid uncontroled effects arising 
from the bubble itself. The Fourier transform of the den- 
sity fluctuation is shown in Fig. [8] at t = 8.5 ps and 
P = bar. While both functionals generate low-g den- 
sity waves in the phonon region (see Figjl]) the ST ap- 
proach does not display any structure, whereas in the OT 
approach one can identify two distinct peaks. The higher 
one is located at q ~ 0.8 A -1 near the maxon region, and 
the lower one is located at q ~ 2.3 A -1 close to the roton 
minimum. 

With these results on the evolution of the IP e-bubble 
in mind, we can state with some confidence that the ST 
description is accurate enough for describing the fate of 
the e-bubble, yielding the appropriate final topology at 
a given pressure, and a more than qualitative picture of 
its evolution. The shape of the cavity, which is respon- 
sible for most electron properties, is essentially the same 
in both ST and OT descriptions. The different way of 
energy release associated with their each kind of elemen- 
tary excitation may yield somewhat diverse evolutions at 
longer times, but it has little relevance for the problems 
at hand. One should keep in mind however, that if the 
actual subject of the study are the elementary excitations 
of the bulk liquid, the use of the OT functional is then 
unavoidable. 

We also want to stress that computing the IP e-bubble 
dynamics in real time for the ST functional has allowed 
us to explicitly check the adiabatic approximation in the 
electron evolution during the time interval relevant for 
the e-bubble "fission" We have computed the overlap 
between the time-evolving electron wave function and the 
instantaneous eigenstate la - , and have found it to be 
equal to unity at all times, meaning that the adiabatic 
approximation holds. As we will discuss later on, this is 
not the case for the 2P e-bubble. 



2. Time-resolved absorption spectrum 

Within the OT functional plus adiabatic approxima- 
tion scheme, we have studied the excitation of IP bub- 
bles by photoabsorption either to the m = component 
(2<r + ), or to the m = ±1 components (l7r - ), arising from 
the splitting of the originally spherical ID state, see Fig. 
|2j In principle, this can be measured in a pump-probe 
experiment by which the e-bubble is excited by two con- 
secutive laser pulses. The delay set between these pulses 
corresponds to the time interval between the excitation 
and the measurement, which is the same as the time de- 
fined in our calculations. The intensity of the absorption 
lines is characterized here by their oscillator strength cal- 
culated in the dipole approximation 3 -^ 
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f ab = 2 ^(E a -E b )\(a\r\b)\ 2 . (7) 

We recall that this oscillator strength fulfills the sum rule 
^2 a f ab = but * s generally not positive-definite. If the 
initial state is not the ground state, a partial sum may 
be greater than unity. 

Starting from the IP electron state la~ (m = 0), the 
two possible photoexcitation transitions are la~ — ^ 2a + 
and la~ — » ±7r _ , see Fig. [2] The specularly asymmetric 
states have a nodal point on the z = symmetry plane, 
implying that they are rather insensitive to the presence 
of helium in that plane. Therefore, the transition energy 
for excitations between two asymmetric states should not 
depend much on whether the bubble has split or not. 
Contrarily, the specularly symmetric states do not have 
such a nodal point, and so are more sensitive to splitting. 
The lowest-lying transition connecting specularly asym- 
metric with specularly symmetric states may thus probe 
the topology of the bubble, since the absorption energy 
for this transition should increase by a sizeable amount 
when the bubble splits. The level structure at the right 
part of Fig. [2] may help understanding these issues. 

The time-resolved absorption energies and oscillator 
strengths of the evolving IP bubble for P = 0.5 and 1 bar 
are presented in Figs. [9]and[l0j As shown by our calcu- 
lations, the IP bubble does not split for P = 0.5 bar, but 
it does for P = 1 bar at t ~ 170 ps. The bubble splitting 
yields a clear signature in the energies and the oscillator 
strengths: the evolution of the transition energies is sim- 
ilar for P = 0.5 and P = 1 bar before the splitting, but 
when the bubble splits at P = 1 bar, the 1<j~ — >• 2<j + 
energy rapidly increases by ~ 70 meV, becoming compa- 
rable to the lcr~ — » 1tt~ energy. This is a consequence of 
the change in the bubble topology, which makes the final 
symmetric and antisymmetric states nearly degenerate. 

A conspicuous pattern also appears in the evolution 
of the oscillator strength. The strength for the specu- 
larly asymmetric transition la~ — >• lir~ remains nearly 
constant at / ~ 0.65, whereas the strength for the 
specularly symmetric transition 1<j~ — >> 2<j + oscillates 
when the bubble is simply connected but falls down to 
/ ~ 0.32 when the bubble splits. This is again a conse- 
quence of the near degeneracy of symmetric and asym- 
metric states in the split-bubble regime. The oscillator 
strength for the antisymmetric transition is a factor of 
two larger than that of the symmetric transition in the 
split-bubble regime because the final state 1tt~ is twofold 
degenerate. 34 

We thus conclude that time-resolved absorption ener- 
gies are of practical interest because they bring rich in- 
formation on the bubble shape and can be determined 
in experiments. This may shed light on the longstand- 
ing question about whether IP e-bubbles under pressure 
do really "fission" into two baby bubbles as our calcu- 
lations indicate, and how the electron wavefunction col- 
lapses into one of them, without violating the quantum 
measurement axiom. 



B. 2P e-bubble dynamics: the breakdown of 
adiabaticity 

Our previous analysis of the dynamics of the IP e- 
bubble has shown that one does not need to use the ac- 
curate OT functional to describe this process. The much 
simpler ST approach already yields a fair description. 
This is particularly useful when we move to the study of 
the 2P e-bubble dynamics. An attempt to simulate this 
evolution has been made within the OT approach and the 
adiabatic approximation!^ This could only be performed 
for a few picoseconds, as it was shown that the adiabatic 
approximation fails at t ~ 7.4 ps. This failure is due to 
the approaching of the 3<j~ (arising from the spherical 
2P level) and 2cr~ (arising from the spherical IF level) 
energy levels. By using the efficient real-time ST scheme, 
we now relax the adiabatic approximation, following the 
evolution of the 2P e-bubble for several hundreds of pi- 
coseconds. We keep refering to this bubble as a "2P e- 
bubble" , but should have in mind that once the adiabatic 
approximation breaks down, the electron is no longer in 
the original eigenstate. Generally, it is in a superposition 
of states that have the same quantum numbers as the 
initial state, meaning that it can be in any superposition 
of a~ states. 

For the 



The 2P bubble evolution is shown in Fig. 11 



first 100 ps, the shape evolution of the 2P bubble is sim- 
ilar to that of the IP bubble, as it expands along the 
symmetry z-axis while its waist shrinks in the perpen- 
dicular plane. From this point on, the bubble oscillates 
back and forth in a kind of four-lobe shapes quite differ- 
ent from those seen in the IP bubble. We attribute these 
conspicuous shape variations to the breaking down of the 
adiabatic approximation as the electron moves from a 
eigenstate to a nontrivial superposition of those compat- 
ible with the symmetries of the system. After evolving 
for ~ 325 ps, the 2P bubble splits into two baby bubbles. 

We present in Fig. 12 'a) the evolution of the instan- 
taneous eigenenergies of the first a~ states. As can be 
seen in panel (b), the 2a~ and 3a~ states nearly meet at 
t ~ 7.4 ps. In agreement with some well-known results 
from basic quantum mechanics J 3 -^ we have found that this 
situation corresponds to an avoided crossing. Panel (c) 
shows the overlap of the evolving electron wavefunction 
with the relevant instantaneous eigenstates. The electron 
is initially in a 3<j~ state (the overlap is unity), but at 
the point of avoided crossing the adiabaticity is lost: the 
electron state is a superposition of the 2a~ (~ 80%) and 
3a~ (~ 20%) states. 

We have also found a time interval between ~ 155 and 
~ 180 ps when the 2P bubble at P = bar has split but 
the emerging baby bubbles do not go away. When this 
happens, the ncr~ and n<r + states should be degenerate. 
This is illustrated in panel (a) of Fig. 12 for the n = 1 
states. Notice from Fig. [IT] that in the 50 ps < t < 100 ps 
interval the bubble is simply connected and the apparent 

'c) is due to the energy 



degeneracy displayed in Fig. 12 



scale. The same thing happens around t ~ 230 ps. 
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IV. SUMMARY 

We have thoroughly studied the dynamical evolution of 
IP and 2P excited electron bubbles in superfluid 4 He at 
zero temperature. To this end, we have resorted to zero- 
and finite-range density functionals, establishing how re- 
liable the former is by comparing its results with those 
obtained with the latter. 

Although the results obtained for the IP bubble evolu- 
tion from these two functionals show some quantitative 
differences, especially for long-time evolutions, they are 
qualitatively equivalent. In particular, both lead to the 
conclusion that IP bubbles "fission" at pressures above 1 
bar. The ST functional result is of particular relevance, 
as it has been obtained by a real-time evolution, with- 
out assuming the adiabaticity of the process. This con- 
firms the previous results obtained using the finite-range 
OT functional and the adiabatic approximation for much 
shorter periods of time than in the present work.^ 

Some indirect experimental evidence indicates a 
change in the de-excitation behavior of the IP e-bubble 
as pressure increases ! 1 2 1 1 6 1 We have explored here the pos- 
sibilities offered by the photoabsorption spectrum of the 
IP e-bubble to disclose whether such a bubble de-excites 
by "fission" or by a more conventional radiative decay, 
and have obtained the signatures that would help distin- 
guish between both decay channels. Although far from 



trivial, a pump-probe experiment may detect a change 
in the absorption spectrum of the IP bubble associated 
with the appearance of the two baby bubble de-excitation 
channel. 

Finally, we have studied the evolution of the 2P e- 
bubble in real-time within the ST functional approach. 
We have dynamically found that the adiabatic approxi- 
mation does not hold at any positive pressure confirm- 
ing the results obtained within the OT plus adiabatic 
approximation approach.^ Negative pressures, as those 
attained in cavitation experiments, have not been stud- 
ied. The physical realization of a relaxed, metastable 2P 
configuration is discarded. 
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FIG. 1: (Color online) Dispersion relation of the elementary 
excitation in bulk liquid 4 He at T = 0. Solid line: the OT 
finite-range functional results. Dashed line: the ST zero-range 
functional results. Dots: the experimental data from Ref. 1221 
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FIG. 2: (Color online) Splitting of electronic levels along with 
a representation of their probability densities once the spher- 
ical symmetry is broken. The states in the spherical configu- 
ration (left) are labeled in the standard nL way. In the axi- 
ally symmetric case (right) the label is nl s z , where s = +( — ) 
stands for symmetric (antisymmetric) states under specular 
reflection. 
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FIG. 3: (Color online) Evolution of the IP e-bubble at P = 
bar. The left-hand side of each panel shows the results for the 
OT finite-range density functional plus adiabatic approxima- 
tion for the electron. The right-hand side part of each panel 
shows the results for the ST zero-range density functional plus 
real-time evolution for the electron. 
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FIG. 4: (Color online) Same as Fig. [3]at P = 0.5 bar. 
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FIG. 5: (Color online) Same as Fig. [3] at P = 1 bar. 
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FIG. 6: (Color online) Same as Fig. [3] at P = 5 bar. 
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FIG. 7: (Color online) Evolution of the extent of the bubble 
along the z axis at P — and 1 bar. The solid line is the OT 
finite-range result, and the dashed line is the ST zero-range 
result. 
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FIG. 8: (Color online) Fourier transform of the density fluc- 
tuation along the z-axis within the region 30 A < z < 70 A 
for the expansion process of the IP e-bubble at P — bar 
and t = 8.5 ps. The solid line corresponds to the OT calcu- 
lation, and the dashed line corresponds to the ST functional 
calculation. The lines have been drawn as a guide to the eye. 
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FIG. 9: (Color online) Time-resolved absorption energies at 
P — 0.5 and 1 bar for the IP e-bubble evolution. The thin 
vertical line at t = 170 ps indicates the time at which the 
bubble splits in the P = 1 bar case. 
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FIG. 10: (Color online) Time-resolved absorption oscillator 
strengths at P — 0.5 and 1 bar for the IP e-bubble evolution. 
The thin vertical line at t — 170 ps indicates the time at which 
the bubble splits at P = 1 bar 
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FIG. 11: (Color online) Evolution of the 2P e-bubble at P = 
bar using real-time dynamics and the ST zero-range func- 
tional. 




FIG. 12: (Color online) (a): Lower-lying instantaneous a~ 
eigenstates together with the lcr + eigenstate of the 2P e- 
bubble at P — bar as a function of time. The thin vertical 
line at t = 325 ps indicates the time at which the bubble 
splits, (b) Enlarged view of the region where the 3a~ and 
2<j~ states repel each other, (c) Overlap of the time-evolving 
electron state onto the 3a~ (solid line) and 2<j~ eigenstates 
(dashed line), \{&(r,t)\na~)\ 2 . The adiabatic approximation 
fails if the value of this overlap varies in time. 



